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The effect of a meridional flow on Parker's interface 
dynamo 
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ABSTRACT 

Parker's interface dynamo is generalized to the case when a homogeneous flow is 
present in the high-diffusivity (upper) layer in the lateral direction (i.e. perpendicular 
to the shear flow in the lower layer) . This is probably a realistic first representation of 
the situation near the bottom of the solar convective zone, as the strongly subadiabatic 
stratification of the tachocline (lower layer in the interface dynamo) imposes a strong 
upper limit on the speed of any meridional flow there. 

Analytic solutions to the eigenvalue problem are presented for the cases of vanish- 
ing diffusivity contrast and infinite diffusivity contrast, respectively. Unlike the trivial 
case of a homogeneous system, the ability of the meridional flow to reverse the propa- 
gation of the dynamo wave is strongly reduced in the interface dynamo. In particular, 
in the limit of high diffusivity contrast relevant to the solar case it is found that a 
meridional flow of realistic amplitude cannot reverse the direction of propagation of 
the dynamo wave. The implications of this result for the solar dynamo problem are 
discussed. 
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1 INTRODUCTION 

While the o peration of the solar dyna mo is still far from 
understood Jweisdll994t iPetrovavfe oO0 1 . it is now generally 
believed that the strong toroidal magnetic field responsible 
for solar activity is generated and stored in the tachocline 
layer. This transitional layer between the differentially ro- 
tating convective zone (CZ) and the rigidly rotating solar 
interior is characterized by a strong radial shear, and is thus 
an ideal candidate for the production of the toroidal field. 

The site and physical nature of the toroidal— >poloidal 
flux conversion (a-effect), needed to close the cycle, is 
much less clear. One popu lar possibility is the inter- 
face dynamo l| Parker| Il993t— hereafter P93lTobiad 1 1996 



Charbonn eau fc MacGregorlll997l . lMarkiel fc Thomas ^ 1998 
Mason et alJl2002T) where a is concentrated in the deepest 



part of the CZ, immediately above th e layer of str ong shear. 

Helioseismic inversions jKosovichevI Il996t 
iBasu fc Antial l200ll) indicate that at low latitudes the 
solar tachocline lies immediately below the adiabatically 
stratified CZ. This implies that the turbulent diffusivity in 
the tachocline is significantly suppressed in comparison to 
the layers immediately above it, resembling the situation 
envisaged in interface dynamo models. An attractive 
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feature of these models is that the sharp diffusivity contrast 
between the two layers gives rise to quite strong toroidal 
magnetic fields, in agreement with the requirements of 
flux emergence calculations. If the sign of a is negative on 
the northern hemisphere, as expected near the bottom of 
the CZ, then the Parker- Yoshimura sign rule predicts an 
equatorward propagating dynamo wave at low latitudes 
and a poleward propagating wave at high latitudes, with 
the two belts departing at the corotation latitude of ~ 35° . 
It is tempting to identify these waves with the two well 
known branches of the exte nded but terfly diagram (e.g. 
[Makarov fc Sivaramanl Il989t cf. also iPetrovav fc Sz akalv 

The presence of a meridional circulation, however, sig- 
nificantly complicates the picture. As a poleward meridional 
flow of amplitude ~ 10 -20 m/s is clearly detected in the shal- 
lower layers of the CZ jHathawavll9 96'l , continuity seems to 
require a counterflow in the deep CZ (though a two-cell pat- 
tern is at present also not excluded). As this velocity ampli- 
tude is comparable to the speed of migration of the activity 
belts, the question arises whether the meridional flow can in- 
validate the Parker- Yoshimura sign rule, reversing the prop- 
agation of a dynamo w ave. Indeed, in flux transport mod- 
els of the solar dynamo ( Wan g et al 1 ll99lHChoudhuri et al. 
Il995lh a deep equatorward meridional flow is responsible for 
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the migration of sunspot-forming latitudes during the solar 
cycle. 

As the effect of a homogeneous flow can obviously be 
described by a Galilean transformation of the solution with 
no flow, it is indeed to be expected that an equatorward 
meridional flow pervading the whole dynamo region near 
the bottom of the CZ (as assumed in flux transport mod- 
els) will reverse the poleward propagation of the dynamo 
wave if its speed is higher than the phase velocity. (In these 
models, a is concentrated near the surface and is positive in 
the northern hemisphere, leading to poleward propagation 
at low latitudes — hence the need for a reversal.) The actual 
situation is, however, more complex, as a meridional flow of 
significant amplitude cannot be expected to penetrate be- 
low the adiabatically stratified CZ. The strong subadiabatic 
stratification of the tachocline represents a serious obsta- 
cle in the way of meridional circulation. The timescale of 
any meridional flow here cannot be shorter than the (tur- 
bulent) heat diffusion timescale, allowing downmoving fluid 
elements to get rid of their strong buoyancy. Heat diffu- 
sivity in the subadiabatic layer is clearly much lower than 
above, so the same must be true for the amplitude of the 
poloidal flow. A more realistic representation of the merid- 
ional flow pattern in the dynamo layer, then, is an interface- 
type model with no meridional flow in the strongly sheared, 
low-diffusivity lower layer, and a homogeneous meridional 
flow in the highly diffusive upper layer. 

This paper addresses the problem of how a meridional 
flow influences the properties of the dynamo wave by con- 
sidering the arguably simplest nontrivial case: the effect of a 
homogeneous meridional flow in the top layer of Parker's in- 
terface dynamo. Section 2 describes the analytic model, Sec- 
tion 3 presents some solutions for important special cases. 
Finally, in Section 4 the implications of our findings are dis- 
cussed for the problems outlined above. 



2 MODEL AND DISPERSION RELATION 

Our model is a generalization of Parker's Cartesian inter- 
face dynamo (P93). The model describes dynamo activity 
occurring in two adjacent layers, say z < and z > 0. This 
Cartesian model can be considered a local approximation 
of the true spherical system, with the x and y coordinates 
corresponding to colatitude and longitude, respectively. The 
magnetic field is written in the form 

B = Be y + V x (Ae y ), 

while the zonal flow is written in terms of a shear rate Q: 

Uy = QZ 

The parameters are homogeneously distributed in each 
layer, and change discontinuously across the interface. In 
the upper part, we assume u x = uo, £1 = 0, a — ceo, V ~ 
while in the lower part, we have £1 = Qo, a = u x = 0, 
T) = r]i . (Here and in what follows, the indices j and J, refer 
to values in the z > and z < sublayers, respectively.) 
The aSl-dynamo equations in the two layers read 



d t A + u d x A = 77V 2 A + ceB, 
d t B + u d x B = r)V 2 B + Qd x A 



with the appropriate boundary conditions ensuring conser- 
vation of magnetic flux applied at the interface (z — 0). 

Following P93, the solution of equations (0 and is 
now sought in the form 



(3) 



Bi = C exp(at — Sz) exp[i(uit + kx — Qz)] 
By = (L + Mz) exp(<rt + sz) exp[i(ujt + kx + qz)] 
Ay — (D + Ez) exp(at ~ Sz) exp[i{ujt + kx — Qz)] 

Ay — J exp(at + sz) exp[i(u}t + kx + qz)] 

where all parameters are real. At the interface z = 0, we 
require 

Ai = Af , dzA t = d z Af, B l =B u = m d z B v 

Nondimensional parameters are now introduced as 



diffusivity contrast: 
dynamo number: 
Reynolds number: 
dimensionless frequency: 
dimensionless growth rate: 



m 2 = mlm, 

N = ce Sl /rfk 3 , 

R = uo/r]^k, (5) 

v = ui/rifk 2 , 

P-l = a/ m k 2 . 



Together with equation this implies that the nondimen- 
sional phase speed of a mode is —v. 

Substituting the ansatz JHJ into equations (4) and (5), 
taking into account the conditions @, and following the 
procedure of P93, we find the complex dispersion relation 

Z 8 ^ 2 - f) 2 + 2Z 6 (^ 2 - l) 2 {n 2 -1-iR) 
+Z 4 [(1 - 3 M 2 + M 4 )(M 2 - 1 - iR) 2 ~ »(1 + 
-zV[(m 2 - l - iR? + »(l + M 2 )0 2 - 1 - iR)N] 



-H 4 N 2 /16 = 0, 



(6) 



(1) 
(2) 



with Z 2 = f3 + i(v + R). By separating equation JHJ into real 
and imaginary parts, one arrives at the dispersion relations 
to be solved for the real quantities (3 and v. While these rela- 
tions can be given in a closed form, they are rather lengthy 
and awkward to handle in the general case. However, the 
dispersion relations simplify significantly for some impor- 
tant special cases. In what follows, we will consider each of 
these cases in turn. 



3 SOLUTIONS FOR SPECIAL CASES 
3.1 Results for the Case /_t = I 

The real and imaginary parts of the dispersion relation © 
now simplify to 

P 2 R 2 + (3N(R + 2v) - vR 2 (R + v) - N 2 /16 = 0, (7) 

(3 2 N + [3R 2 (R + 2v) - vN{ R + v) = 0. (8) 

Multiplying equation Q by N and equation @ by R 2 , and 
taking the difference we have 

j3{R+2v){R i + N 2 ) - 7V 3 /i6 = 0. (9) 
From this, for the frequency v we have the explicit result 
iV 3 R 



32f3(R 4 + N 2 ) 2 



(10) 
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Figure 1. Dynamo wave period P = 2tt/\v\ (solid), phase speed 
—v (dotted), and growth rate /3 — 1 (dashed), as functions of the 
flow speed R in the analytic model with uniform diffusivity fj, = 1, 
with N = -100. 



while for the growth rate /3 
P 2 : 



1 we a quadratic equation in 



N 6 - 64jV 4 (16/3 4 + 3R 2 (3 2 ) 
-256(47? 8 /? 4 + R 10 P 2 ) = 0. 



UN 2 (32R 4 (3 4 + 7R 6 /3 2 ) 
(11) 

In Figure the wave period P = 2n/\v\, the phase 
speed —v, and the growth rate 13—1, computed from rela- 
tions ([113 and fTO . are plotted for N = -100. The critical 
dynamo number for R = is iVcr = ±32, so this corresponds 
to a moderately supercritical case. It is apparent that for 
high negative values of R the phase speed becomes slightly 
negative, indicating that the dynamo wave is reversed. Nev- 
ertheless, the minimal flow speed needed for this, R ~ —5, 
is significantly higher in modulus than the phase speed for 
R — (— v ~ 1), in contrast to the case of a uniform flow. 
This confirms that with the nonuniform flow considered here 
the properties of the dynamo differ significantly from the 
case of a uniform flow that simply advects the dynamo wave. 
In particular, the Parker-Yoshimura sign rule 



sgn v = sgn N, 



(12) 



that determines the sense of dynamo wave propagation when 
R = 0, is not as easy to violate as in the case of a uniform 
flow. This is further borne out in Figure^ where the wedge- 
shaped region marked "reverse wave" denotes that portion 
of the N-R plane in which the sign rule is violated. For a 
uniform flow, this region should clearly extend down to the 
dotted line indicating the negative of the phase speed in the 
R = case. 

It is apparent that the growth rate is maximal for R — 
0, and otherwise independent of the direction of the flow 
(as only even powers of R appear in equation llll ). This 
symmetric behaviour of the (3(R) curve is, however, specific 
to the case fj. = 1. Indeed, expanding the dispersion relation 
in terms of m = fi 2 — 1 and R around the case m = R = 0, 
after considerable algebra one finds 



dR 



-3mN 
1024^ 



Figure 2. Minimum antiparallel flow speed necessary to reverse 
the dynamo wave (solid) in the case of uniform diffusivity in the 
analytical model. Dashed line indicates the critical mode sepa- 
rating the sub- and supercritical regimes. The dotted line shows 
the negative of the phase speed — v for R = 0. In the case of a 
uniform flow the solid and dotted lines would coincide. 



(3q — 1 being the growth rate for m = R = 0. From this it 
is clear that for iV < 0, in models with /1 2 < 1, the growth 
rate is increased by a finite negative R, while with /1 2 > 1 



the reverse is true. 



3.2 Results for the Case (i<l 

We now turn to the case more relevant to the Sun, when the 
shear layer is characterized by a lower value of the diffusivity 
than the layer above. Again following P93, we consider the 
limit fj, — > while N = fi 2 N remains finite, and restrict our 
attention to unstable modes with fi 2 <C (3 — 1. The real and 
imaginary parts of the dispersion relation © can then be 
written as 



(y + R) 4 - 2R(y + R) 3 - {v + R) 2 (6/3 2 -6/3 + 1 
+(u + R)[2/3R(3f3 - 2) + N{fi - §)] + f3 2 (/3 - l) 2 
-f3R{f3R+\N)-±N 2 =0, 



R 2 ) 



(13) 



{{v+R) 2 -(3(P-l)-(v+R)R}[N-8(v+R)(f3- 1 :)+4[3R] = 0.(14) 



The latter relation implies that either 
u=^{-R±[R 2 + 4((3~l)/3} 1/2 }, 



(15) 



N = 8(v + R)(l3-±)-4l3R 



(16) 

Substituting either of the last two relations into 1)13^ yields 
the other, so by combining them we obtain a relation that 
gives the growth rate implicitly 

N = 4{R ± [R 2 + 4{J3 - l)/3] 1/2 }((3 - |) - 4(3R. (17) 

Upon solving this for f3, v can be calculated from equation 
1161 . It is clear that for iV > it is the + sign, while for 
N < it is the — sign that corresponds to growing modes 
in the case R — 0. The other sign will not lead to growing 
modes for any value of R. 
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Figure 3. Dynamo wave period P = 2tt/\v\ (solid), phase speed 
— v (dotted), and growth rate (3 — 1 (dashed), as functions of 
the flow speed R in the analytic model in the limit of very high 
diffusivity contrast u 2 < /3 - 1, with TV = -12\/2. 

The dependence of the wave period P = 2it/\v\, the 
phase speed — v, and the growth rate {3 — 1 on 7? is plotted 
in Figure [3] for the growing mode (— sign), in a moderately 
supercritical case with a negative dynamo number. 

Inspection of Figure[3]and, more importantly, the anal- 
ysis of equations 1151 - dlT| also indicates that — v > for all 
values of R for growing modes; that is, the sign rule 11211 
seems to hold unconditionally in the limit considered here. 
Indeed, while a parallel flow of relatively moderate speed 
-Rmax (~ 4 in the case considered) can quench the dynamo 
completely, an antiparallel flow (R < 0) can apparently nei- 
ther quench the dynamo wave, nor turn it around, for any 
finite speed. However, a caveat is in order when applying 
these findings to the case of a small but finite (i. As the re- 
sults were derived assuming /i 2 -C /3 — 1, their applicability 
is limited to a limited range in R, the limits of which are 
determined by the condition f3(R) - 1 > /i 2 . However, from 
equation 1171 in the case of an antiparallel flow one finds 
the condition R > R CI = iV/4/i 2 (for N < 0), the order of 
magnitude of R cv being at least 100, much higher than any 
realistic estimate of R in the solar case. Thus we can safely 
conclude that in the case of high diffusivity contrast and 
meridional flow limited to the upper regime ( CZ), a merid- 
ional flow cannot reverse the direction of propagation of the 
dynamo wave, unless its amplitude is unrealistically high. 



a maximum, the growth rate tends to zero as uo — > — oc. 
Contrary to intuition, an antiparallel flow can neither sup- 
press nor revert the direction of propagation of the dynamo 
wave. 

This conclusion is appa rently at odd s with 
flux transport mo d els llChoudhuri et all Il99l 
iDikpati fc Charbonneaul Il999l) where an equatorward 
flow near the bottom of the CZ can turn the direction of 
propagation of the dynamo wave towards the equator at 
low latitudes. There is, however, no real contradiction here, 
as in those models the flow is assumed to penetrate quite 
deep into the layer of strong shear. In terms of the simple 
Cartesian model studied here, a closer counterpart of the 
flux transport models would be a case with uo =const. 
throughout the region: then an antiparallel flow exceeding 
the phase speed can trivially turn around a dynamo wave. 
In Section 1 we argued that it is more realistic to assume 
that the meridional flow is limited to the adiabatic upper 
layer. In fact , the more recent flux transport model of 
IDikpati et alJ i2004h admits that the circulation must be 
limited to the adiabatic layer. Therefore, in order to make 
the model work, they need to assume that a significant 
fraction of the tachocline overlaps the adiabatic SCZ. This 
may indeed hold at high latitudes, but at low latitudes, 
where a meridional transport of toroidal fields is most 
neede d, it does not seem to be supported by helioseismic 
data feasu fc Antiall200ll) . 

These arguments indicate that our simplified models 
may have more relevance to the solar dynamo than many 
complex nonlinear spherical models where the assumed ge- 
ometrical distribution of the flows and transport coefficient 
does not reflect the situation in the solar interior. Never- 
theless, owing to the simplified 2D Cartesian geometry and 
the kinematical nature of our models our conclusions should 
be taken with proper reservation. In this context it may 
be mentioned that that details of the physical structure, 
such as the precise form of the rotation law, were shown to 
hav e a profound effect on the behaviour of d ynamo mod- 
els (iMoss fc Brooke! 120001 . iPhillips et al.ll2002l) . Extensions 
of this work to spherical geometry and to the 3D nonlinear 
domain are clearly needed to tell whether the present results 
remain valid in more realistic situations. 
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4 CONCLUSION 

We have studied by analytical methods how a meridional 
flow limited to the upper (i.e. unsheared) layer affects kine- 
matic interface dynamos. We found that the growth rate and 
period of the dynamo wave have a nontrivial dependence on 
the flow amplitude. In the case of strongly reduced diffusiv- 
ity in the lower layer, relevant for the Sun, a flow parallel 
to the direction of propagation of the dynamo wave reduces 
the growth rate until, at a finite critical flow speed, it com- 
pletely suppresses dynamo action. An antiparallel flow, in 
contrast, first increases the growth rate; then, after reaching 
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